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We reexamine whether ^"'"(4430) could be a D[ — D* or Di ~ D* molecular state after consid- 
ering both the pion and a meson exchange potentials and introducing the form factor to take into 
account the structure effect of the interaction vertex. Our numerical analysis with Matlab package 
MATSLISE indicates the contribution from the sigma meson exchange is small for the D'l — D* 
system and significant for the Di — D* system. The S-wave Di — D* molecular state with only 
= 0~ and Di — D* molecular states with — 0~, 1~, 2~ may exist with reasonable parameters. 
One should investigate whether the broad width of D'^ disfavors the possible formation of molecular 
states in the future. The bottom analog Zb of Z"'"(4430) has a larger binding energy, which may be 
searched at Tevatron and LHC. Experimental measurement of the quantum number of ^"''(4430) 
bJ[), may help uncover its underlying structure. 
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^ ' I. INTRODUCTION 

I ■ 

Recently Belle collaboration observed one very exotic resonance Z+(4430) in the tt+i/;' invariant mass spectrum in 
(— I ■ the exclusive B — > Kn^^' decays Its mass and width are 

M = 4433 ± 4(stat) ± l(syst) MeV 

(N ■ 

J> ' and 

0\ . r = 44j:J^(stat)i5?(syst) MeV. 

CN 

This resonance appears as an excellent candidate of either the multiquark state or the molecular state and has 
stimulated many theoretical investigations 0, i, i, d, MM, S, i, 0, [HI,!!!, [H Q d, [H, El ■ A concise review of 
[ the theoretical status of Z+ (4430) can be found in Ref. [4I . Thereafter Ding studied Z+ (4430) using the approach of 
00 , effective Lagrangian and predicted a B*Bi bound state with mass 11048.6 MeV [l^l- Later Braaten and Lu studied 
O ■ the hne shapes of Z+(4430) 

^ , In our previous work 0, we explored whether Z+(4430) could be a loosely bound S-wave molecular state of D* 
r ' and D'^ (or Di) with = 0~, 1~, 2~. We considered the one-pion exchange potential only and did not introduce the 
rS form factor in the scattering matrix elements in the derivation of the potential. Instead of solving the Schrodinger 
' equation numerically, we employed some simple trial wave functions and the variation method to study whether a 
■ " " ' shallow bound state exists. We found that the interaction from the one pion exchange alone is not strong enough 
to bind the pair of charmed mesons with realistic coupling constants. Other dynamics is necessary if Z+(4430) is a 
molecular state. 

The one pion exchange potential alone does not bind neutron and proton to form the deuteron in nuclear physics 
either. The strong attractive force in the intermediate range is required in order to bind the deuteron, which is modelled 
by the sigma meson exchange. Meanwhile hadrons are not point-like objects. The cutoff should be introduced to 
describe the structure effect of the interaction vertex. With the above considerations, we have reexamined the DD* 
molecule picture for X(3872) after taking into account both the pion and sigrna meson exchange. It turns out that the 
sigma meson exchange potential is repulsive and numerically important [l8|. One may wonder whether the similar 
mechanism plays a role in the case of Z+(4430). Therefore we will make a comprehensive study whether Z+(4430) is 
the molecular state. 
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This work is organized as follows. We present the effective Lagrangian relevant to calculate the tt and a exchange 
potentials for Z+(4430) in Section [TTl We present the detail of the derivation of the potential in Section Hill We 
discuss the behavior of the potential in Section Hvl We make numerical analysis in Section |Vl We discuss the bottom 
analog of Z+(4430) in Section IVTl We discuss the cutoff dependence in Section IVlII The last section is the conclusion. 



II. THE EFFECTIVE LAGRANGIAN AND COUPLING CONSTANTS 



In order to derive the tt and a exchange potentials, we collect the relevant effective Lagrangian in this section. The 
Lag rangian for the interaction of tt and charmed mesons is constructed in chiral symmetry and heavy quark limit 

Ct, = igTT[Hb4-balbHa] + ig'TT[Sb4balbSa] + ig"TT[T^ib4balbT^] + ihTT:[Sb4balbHa] 

+l!^i:i[T^{D^4)balM + t^TT[T^^mf.)bal5Ha] + h.C. (1) 

where the fields H, S and T are defined in terms of the (0^, (0+, 1+), 2+) doublets respectively 

Ha = ^[P:'l,-Pal,], (2) 
Sa = ^[Kl,l^-Poa], (3) 

Here the axial vector field A'^i^ is defined as 



Kb = lii^d^i - ^d''^^)ab = j^d'Mab 



with ^ = ex^piiM/U), L ^ 132 MeV and 



^ + ^ TT+ K+ 



M = \ . -^ + -^ A- ]. (5) 

K- K° 

In Eq. JT]) , the coupling constants were estimated in the quark model fl9\ , 

9 = gA, g' = 5^/3, g" = gA (6) 

with gA = 0.75. A different set of coupling constants can be found in Ref. [2l|. With our notation, we get the 
following relations 

g = g'A, g' = -3a. h = ^Ga (7) 

with Ga ~ 1 and q'a = 0.6 12111 . In fact, the coupling constant gA was studied in many theoretical approaches such 
as QCD sum rules |22l. [2^ |24|. [25| . In this work we use the experimental value gA — 0.59 ± 0.07 ± 0.01 extracted 
from the width of D* '2u\. With the available experimental information, Casalbuoni and collaborators extracted 
h = 0.56 ± 0.28 and h' = {hi + h2)/k^ = 0.55 GeV"^ |20|]. The signs of g, g' g" are not determined affhough their 
absolute values are known. 

The interaction Lagrangian related to the a meson can be written as 



Tv[S{d^a)-i^H] + -fi:i[T^'{d^a)H] + h.c. (8) 



= g,i:i[HaH]+g'„TY[S(jS]^g,i:T[T^aT,,] 

Jtt Jit 

In order to estimate the values of the coupling constants, we compare the above Lagrangian with that in Ref. [2lj 
and get 

5. = -^5., 9'. = ~:^^g., h. = ^^9'A, (9) 
where gj^ = 3.73. When performing numerical analysis, we take g'^ = g^ and h'^ — approximately. 
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III. DERIVATION OF THE PION AND SIGMA EXCHANGE POTENTIAL 

If Z+(4430) is a molecular state of D[ - D* or Di - D* , the flavor wavefunction of Z+(4430) reads 0, 



|Z+(4430)) = -^\\d[^D*+) + \D*"d'+) 
v2 L 



(10) 



|Z+(4430)) 



V2 



\D°D* 



\D*^Dp 



(11) 



Here D[ and Di with quantum number = 1+ belong to (0+, 1+) and (l+,2+) doublet respectively in the heavy 
quark limit. 

To derive the effective potentials, we follow the same procedure as in Ref. [2]. Firstly we write out the elastic 
scattering amplitudes of the direct process A{B) A{B) and crossed channel A{B) B{A), where A and B denote 

and D*. Secondly, we impose the constraint that initial states and final states should have the same angular 
momentum. Thirdly, we average the potentials obtained with Breit approximation in the momentum space. Finally 
we perform Fourier transformation to derive the potentials in the coordinate space. For the scattering between D* 
and D'l (Di), both the pion and sigma meson exchange are allowed in both direct and crossed processes. 

We introduce the form factor (FF) in every interaction vertex to compensate the off-shell effects of the exchanged 
mesons when writing out the scattering amplitude, which differs from Ref. One adopts the dipole type FF (27l.l28j 



F{q) = 



A2- 



A2 



(12) 



The phenomenological parameter A is near 1 GeV. q denotes the four- momentum of the exchange meson. It is observed 
that as ^ it becomes a constant and if A ^ M, it turns to be unity. In the case, as the distance is infinitely 
large, the vertex looks like a perfect point, so the form factor is simply 1 or a constant. Whereas, as — > oo, the 
form factor approaches to zero, namely, in this situation, the distance becomes very small, the inner structure (quark, 
gluon degrees of freedom) would manifest itself and the whole picture of hadron interaction is no longer valid, so the 
form factor is zero which cuts off the end effects [l^. 

For the direct scattering channel, go is a small value. In the heavy quark limit we approximately take go — 0. Thus 
it is reasonable to approximate q^ as — q^. However, for the crossed diagram, we can not ignore the contribution from 
go due to go = — Mo' ~ 410 MeV, which is about three times of the pion mass. The principal integration is a 
good approach to solve the problem when go is larger than the mass of exchanging meson. 

Since we only consider S-wave bound states between Di(Di) and D*, there are five independent parts related to 
the potentials of D^{Di) — D* system in the momentum space. We use the following definitions to denote them after 
performing Fourier transformation. 



1 



-I- 



9 9 

g^ — m-^ 

2\2 



(q^) 



9 9 

q- — 



9 9 

g^ — TO^ 



>o(A, m, r), 
Yi{A,m,r), 
Y2{A,m,r), 

Y4{A, m, r). 



(13) 
(14) 
(15) 
(16) 
(17) 
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where m denotes the mass of the tt or cr meson. Their exphcit expressions are 

.,'2 

e 



with 



yo(A, mo-,r) = 


47rr 


Yi(A,TO^,r) = 


2 

47rr 


F2(A, m^,r) 


2 

go 

47rr 


Y3(A,TO^,r) = 


47rr 




[e 

47rr 



^'')-^e-^^-^/'e(A)-rfC(A), (18) 
_ ^_A.^ _^ ^e-^-- + m2^4^(A) + A2,y6^(A), (19) 
cos(H - e-"1 + l^^e-"'- + g2^4^(a) + go2r;'^C(«), (20) 

4 2 4 

cos(H + ^e-"'- + li^e-"'- + - 2Ai'a')^(a) + a\^C(a), (21) 

47rr OTra 

87ra 



^^^^ " 192^^ 64^ "^64^' ^^^^ 



where = V 1o ~ ^1 ' = \/ — gg, ry = y^A^ — m^, ?y' = ^ A^ — and a = \J h? — q^. With the above 
functions, the potentials for the different cases are coUected in Table HI and HIl 





VDir(r) 




Vcro 


(r) 






TT exchange 


a exchange 


TT exchange 




fj exchange 


J= 




^Fi(A, m^,r) 




-^F2(A,m^ 


r) 


|yf-F4(A, m<^,r) 


J= 


=1 




g^g'„Yo{A,m^,r) 


-^F2(A,m^ 


r) 


i^F4(A, m,,r) 


J= 


=2 


-ff;Fi(A,m.,r) 


gag'„Yo{A,ma,r) 


-jjjY2{A,m^ 

J TT 


r) ^ 


-k^Y4,{A,m^,r) 



TABLE I: The potentials with the D[-D* molecular assumption for ^"''(4430). 



J=0 
J=2 



VMr) 

TT exchange a exchange 

-fff;Fi(A,m,,r) 3„ff;'Fo(A,m„,r) 

-fff^Fi (A, m^,r-) g„g'J,Yo{A,m„,r) 
|f|^Yi(A,m^,r) g^g'iYo{A,m„,r) 



Vcro{r) 

TT exchange o exchange 

^y3(A,m^,r) i(^F4(A, m,,r) 

-^^Yi{A,m^,r) ±^J^Y4{A,m^,r) 

^Ys(A,m^,r) (^^y4(A, m., r) 



TABLE II: The effective potentials with the assumption of Z^(4430) being a Di-D* molecule. 



Assuming Z+(4430) to be D[D* or DiD* molecule state, the total potential is 

^Total (r) = Vuir (r) + Vcro (0 , (25) 

where the sign between Vuirir) and Vcmir) is determined by the flavor wavefunction of Z+(4430) in Eq. (fTO|) or (jlip. 
^Dir('') and Vcro('") correspond to the potentials from the direct and crossed diagram respectively. From Tables U 
and im we make two interesting observations: (l)Vcro(?') does not depend on the sign of coupling constant; (2)For the 
D[D* systems, the sigma exchange potentials for the direct diagram and the pion exchange potentials for the crossed 
diagram are the same for the three cases with J = 0, 1, 2! 



IV. THE SHAPE OF THE PION AND SIGMA EXCHANGE POTENTIAL 



In this section we study the variation of the pion and sigma meson exchange potentials with the coupling constants, 
which are given in the Section |TT1 We also need the following input parameters rriD* = 2008.35 MeV, m^i' = 2427 
MeV, = 2422.85 MeV; A = 132 MeV, = 135.0 MeV, m„ = 600 MeV 0. 
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A. Single pion exchange potential 



For the S-wave D[D* system, we take several typical values of the coupling constants g ■ g' = ±0.1, ±0.5 and 
h — 0.56, 0.84. Meanwhile we take the cutoff A = 1 GeV. We illustrate the dependence of single pion exchange 
potential on these typical values in Fig. [1] We also plot the potential of the DiD* molecular state with several typical 
coupling constants [g ■ g" , /i'] = [±0.2, 0.55 GeV"!], [±0.6, 0.55 GeVi] in Fig. H 
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FIG. 1: The single pion exchange potential for the D'-^D* molecule. In diagrams (a), (c) and (e), the solid, thick solid, 
dotted and dashed lines correspond to the potentials with parameters \g ■ g' , h] = [0.1, 0.56], [0.5, 0.56], [0.1, 0.84], [0.5, 0.84] 
respectively. In (b), (d) and (f), the solid , thick solid, dotted and dashed lines correspond to the potentials with parameters 
[g ■ g', /i] = [-0.1,0.56], [-0.5, 0.56], [-0.1, 0.84], [-0.5, 0.84] respectively. Here A = 1 GeV. 




r [GeV-'l r [GeV"'] ' [GeV-'] 

(d) (e) (f) 

FIG. 2: The single pion exchange potential for the DiD* molecule. In (a), (c), (e), the solid and dotted lines denote the 
potentials with parameters [g ■ g" , h']=[0.2, 0.55 GeV~^], [0.6, 0.55 GeV~^] respectively. In (b), (d), (f), the solid and dotted 
line denote the potentials with parameters [g ■ g" , h'] — [-Q.2, 0.55 GeV~^], [-0.6, 0.55 GeV~^] respectively. Here A = 1 GeV. 



From Fig. [T] and [21 we notice that (1) the variation oi g ■ g' or g ■ g" does not result in the big change of the 
potential when r is larger than 6 GeV~^; (2) the potential is sensitive to the value of h, which indicates the single pion 
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exchange in the crossed diagram plays an important role to bind the D^D* compared with the contribution of direct 
diagram; (3) for the S-wave DiD* system with J = 0, its potential is repulsive with [g ■ g" , h'] — [-0.6, 0.55 GeV^^]. 
The potentials of the J = 1 DiD* system with [g ■ g" , h'] = [-0.2, 0.55 GeV-^], [-0.6, 0.55 GeY'^] and the J = 2 DiD* 
system with [g ■ g", h'] = [0.2, 0.55 GeV~^], [0.6, 0.55 GeV~^] are also repulsive, which are shown in Fig. [5] (b), (d) 
and (e). In the range r < 6 GeV^^, the potential of the Di — D* system is sensitive to the coupling constants. These 
conclusions were obtained when taking A = 1 GeV. 



B. The potential via both pion and sigma exchanges 



We further investigate the potential of the D'^D* system by adding sigma exchange contribution. The typical values 
of coupling constants are given in the captions of Figs. [3]l5l In these figures, we compare the single pion exchange 
potential with the total potential with different coupling constants. 




(a-1) (a-2) (a-3) (a-4) 

FIG. 3: For the D'^D* molecular state with J — 0, we compare the single pion exchange potential with the total potentials 
containing the sigma exchange. The thick solid lines denote the single pion exchange potentials. The solid lines, thin solid line, 
dashed line and dotted line correspond to the potentials with parameters [gag'^^, hc,] = [0.58, 0.4], [0.58, 0.8], [-0.58, 0.4], [-0.58, 
0.8] respectively, (a-1), (a-2), (a-3) and (a-4) respectively correspond to [g ■ g' — 0.1, h = 0.56], [g ■ g' — 0.5, h — 0.56], [g ■ g' 
= 0.1, h = 0.84] and [g ■ g' = 0.5, h = 0.84]. Here A = 1 GeV. 
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FIG. 4: For the D'lD* molecular state with J = 1, we compare the single pion exchange potential with the total potentials. 
The thick solid lines denote the single pion exchange potentials. The solid lines, thin solid line, dashed line and dotted line 
correspond to the potentials with parameters [gag'a-, ha] — [0.58, 0.4], [0.58, 0.8], [-0.58, 0.4], [-0.58, 0.8] respectively, (a-1), (a-2), 
(a-3) and (a-4) with [g ■ g' = 0.1, h = 0.56], [g ■ g' ^ 0.5, h = 0.56], [g ■ g' ^ 0.1, h = 0.84] and [g ■ g' = 0.5, h = 0.84]. Here 
A = 1 GeV. 




FIG. 5: For the D'lD* molecular state with J = 2, we compare the single pion exchange potential with the total potentials. 
The thick solid lines denote the single pion exchange potentials. The solid lines, thin solid line, dashed line and dotted line 
correspond to the potentials with parameters [gag'a, ftCT]=[0.58,0.4], [0.58,0.8], [-0.58,0.4], [-0.58,0.8] respectively, (a-1), (a-2), 
(a-3) and (a-4) with [g ■ g' = 0.1, h = 0.56], [g ■ g' = 0.5, h = 0.56], [g ■ g' = 0.1, h = 0.84] and [g ■ g' = 0.5, h = 0.84]. Here 
A = 1 GeV. 
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From Figs. [3][5l one finds that adding sigma exchange contribution does not result in dramatic change of the total 
potential VTotai('') of the S-wave D[D* and DiD* molecule system especially when r is larger than 6 GeV~^. 

V. NUMERICAL RESULTS 

Different from the analysis in Ref. in this work we solve the the Schrodinger equation numerically with the 
help of MATSLISE package, which is a graphical Matlab software package for the numerical study of regular Sturm- 
Liouville problems, one-dimensional Schrodinger equations and radial Schrodinger equations with a distorted coulomb 
potential. It allows the fast and accurate computation of the eigenvalues and the visualization of the corresponding 
eigenfunctions ,30,] . 



A. S-wave D'l — D* system 



In Figs. [6] (a) and [7] (a) we show the radial wave function R{r) and function x(r) — rR{r) for the D'^ — D* 
system with J — Q. We list the numerical results with different typical values of coupling constants in Table Hill for 
the D'l — D* system. Here rrms is the root-mean-square radius and r^ax denotes the radius corresponding to the 
maximum of the wave function x(^) of D'^ — D* system. E{h) denotes the binding energy with the corresponding 
cutoff. For example, the notation -6.0(1.5) denotes the binding energy is 6.0 MeV at the cutoff A = 1.5 GeV. 

From the numerical results listed in Table Hill one concludes that the existence of S-wave D'l — D* bound state with 
= 0~, 1~,2^ is possible. With appropriate parameters, one can get a molecular state consistent with Z+(4430). 
Throughout our study, we have ignored the width of heavy mesons. However, the broad width of D'^ around F ^384 
MeV [29| may be a obstacle for the formation of the molecular state, which deserves further study. 

By comparing the results with different sets of parameters, one finds that the sigma exchange interaction induces 
very small effects on the binding energy. Of the parameters g ■ g' and h, the binding energy is sensitive to h, which 
indicates that the crossing diagram from one pion exchange plays an important role in binding D[D*. Numerically, 
large h, small g ■ g' for J = 0, 1 and big g ■ g' for J — 2 are helpful to form a bound state. These observations are 
consistent with the conclusions by analyzing the dependence of the potentials on different coupling constants in the 
previous section. 

In Table HTTl we only give results corresponding to two cutoffs A = 0.7 GeV and A = 1.5 GeV. By comparing the 
binding energies with different A when h = 0.84, one notes that E becomes larger with a smaller A. We come back 
to this point later when discussing the cutoff dependence of the binding energy E. 
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FIG. 6: The radial wave function R{r) for the molecular state with J = of the D'lD* and DiD* system respectively. The 
corresponding parameters are [g ■ g' ~ 0.1, h = 0.84, A=1.5 GeV] and [g ■ g" = 0.2, h' = 0.55 GeV~^, A=2.9 GeV] respectively. 
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0.4 


-22.3(0.7) /-4.8(1.5) 


1.5/2.2 


1.3/1.6 


-23.3(0.7)7-5.2(1.5) I.572.I 1.371.5 


-25.3(0.7)7-6.0(1.5) I.472.O I.271.5 


0.5 


0.84 


n 8 
u.o 




1 Pi /9 9 
1.0/ Z.Z 


1 /I 

1.0/ i.D 


-23.3(0.7)7-5.1(1.5) 1.572.1 1.371.5 


-25.3(0.7)7-6.0(1.5) I.472.O I.271.5 


-0.58 


0.4 


-22.3(0.7)7-5.1(1.5) 


1.572.2 


1.3/1.6 


-23.3(0.7)7-5.4(1.5) 1.572.1 1.371.5 


-25.3(0.7)7-6.3(1.5) I.472.O I.271.5 






0.8 


-22.4(0. 7) /-5. 0(1. 5) 


1.572.1 


1.3/1.5 


-23.3(0.7)7-5.4(1.5) I.572.I 1.371.5 


-25.3(0.7)7-6.3(1.5) I.472.O I.271.5 


n 1 

-U. 1 


u.ou 






Q o/n 7^ 


9 ^ 


1 4 


-3.1(0.7) 


2.3 


1.4 


-2.9(0.7) 


2.4 1.4 






0.58 


0.4 


-3.2(0.7) 


2.3 


1.4 


-3.1(0.7) 


2.4 


1.4 


-2.9(0.7) 


2.4 1.5 


-0.1 


0.56 


n 8 

u.o 
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1 A 


-3.1(0.7) 


2.4 


1.4 


-2.9(0.7) 


2.4 1.5 


-0.58 


0.4 


-3.2(0.7) 


2.3 


1.4 


-3.1(0.7) 


2.4 


1.4 


-2.9(0.7) 


2.4 1.5 






0.8 


-3.2(0.7) 


2.3 


1.4 


-3.1(0.7) 


2.4 


1.4 


-2.9(0.7) 


2.4 1.5 


n 1 

-U. i 


U.04I: 






9/1 7(^0 7^ / ^ Q('1 


1/1/9 
i.^/ z.u 


1 9 /I Pi 
i.z/ 1.0 


-24.5(0.7)7-5.8(1.5) 1.572.0 I.271.5 


-24.1(0.7)7-5.6(1.5) 1.572.1 1.371.5 






0.58 


n A 


94 7('n 7V ^ 7('1 


± .4:7 


1 9/1 


-24.5(0.7)7-5.6(1.5) 1.572.1 I.271.5 


-24.1(0.7)7-5.5(1.5) 1.572.1 1.371.5 


-0.1 


0.84 


0.8 


-24. 7(0. 7) /-6. 5(1. 5) 


1.471.9 


1.271.5 


-24.5(0.7)7-6.0(1.5) 1.572.0 I.271.5 


-24.1(0.7)7-5.8(1.5) 1.572.1 1.371.5 


-0.58 


0.4 


-24.7(0.7)/-D.0(1.5) 


1.472.0 


1.271.5 


-24.5(0.7)7-5.9(1.5) 1.572.0 I.271.5 


-24.1(0.7)7-5.7(1.5) 1.572.1 1.371.5 






0.8 


-24.7(0.7)7-6.5(1.5) 


1.471.9 


1.271.5 


-24.5(0.7)7-6.0(1.5) 1.572.0 I.271.5 


-24.1(0.7)7-5.8(1.5) 1.572.1 1.371.5 


-0.5 


0.56 




- 


-4.2(0.7) 


2.1 


1.4 


-3.6(0.7) 


2.2 


1.4 


-2.4(0.7) 


2.6 1.5 






0.58 


0.4 


-4.2(0.7) 


2.1 


1.4 


-3.6(0.7) 


2.2 


1.4 


-2.4(0.7) 


2.6 1.5 


-0.5 


0.56 


0.8 


-4.2(0.7)/-l.l(1.5) 


2.173.3 


1.471.4 


-3.6(0.7) 


2.2 


1.4 


-2.4(0.7) 


2.6 1.5 


-0.58 


n A 
u.^ 




2.1710.6 1.471.7 


-3.6(0.7) 


2.2 


1.4 


-2.4(0.7) 


2.6 1.5 






0.8 


-4.2(0.7)7-1.6(1.5) 


2.172.8 


1.471.4 


-3.6(0.7) 


2.2 


1.4 


-2.4(0.7) 


2.6 1.5 


-0.5 


0.84 






-26.4(0.7)7-6.8(1.5) 


1.472.0 


1.271.5 


-25.3(0.7)7-6.2(1.5) 


1.472.0 1.271.5 


-23.3(0.7)7-7.6(1.5) 1.572.1 1.371.5 






0.58 


0.4 


-26.4(0.7)7-6.8(1.5) 


1.471.9 


1.271.5 


-25.3(0.7)7-6.0(1.5) I.472.O I.271.5 


-23.3(0.7)7-5.2(1.5) I.572.I 1.371.5 


-0.5 


0.84 


0.8 


-26.4(0.7)7-8.0(1.5) 


1.471.8 


1.271.4 


-25.3(0.7)7-6.2(1.5) I.472.O I.271.5 


-23.3(0.7)7-5.3(1.5) I.572.I 1.371.5 


-0.58 


0.4 


-26.4(0.7)7-7.2(1.5) 


1.471.9 


1.271.5 


-25.3(0.7)7-6.4(1.5) I.472.O I.271.5 


-23.3(0.7)7-5.5(1.5) I.572.I 1.371.5 






0.8 


-26.4(0.7)7-8.6(1.5) 


1.471.8 


1.271.4 


-25.3(0.7)7-6.6(1.5) I.471.9 I.271.5 


-23.3(0.7)7-5.6(1.5) I.572.I 1.371.5 



TABLE III: The numerical results of D[ — D* system with = 0~, 1~, 2~. Here the bing energy E, cutoff A, the root-mean- 
squaro radius rrms and Vmax of the D'^ — D* system are in unit of MeV, GeV, fm and fm respectively. E{K) denotes the binding 
energy with the corresponding cutoff. 
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FIG. 7: The function ~ rR{r) for the J=0 molecular state of the D'lD* and D\D* system, (c) shows the short range 
behavior for the D\D* case. The corresponding parameters for the two systems are [g ■ g' = 0.1, h — 0.84, A=1.5 GeV] and 
[g ■ g" = 0.2, h' = 0.55 GeV-\ A=2.9 GeV] respectively. 



B. S-wave D\ — D* system 



In Tabic IIV( we present the numerical results for the case of S-wave Di — D* system. Unfortunately one fails 
to find solutions with negative binding energy for J — 1,2 using the parameters in that table, which indicates that 
there probably does not exist the S-wave Di — D* molecule with — 1^, 2~ through the pion and sigma exchange 
interactions. However the S-wave Di — D* molecular state with = 0^ probably exists. Thus we only present the 
results for the Di — D* system with J ~ 0. 

When taking cutoff A = 2.9 GeV, one gets E = -1.3 MeV with parameters [g ■ g",h',g^ ■ .g^', /i^] = [0.2,0.55 
GeV""'^, 0.58, 0.2]. Varying parameters g ■ g", g^ ■ g'l, h'^ in the reasonable range results in large change for the binding 
energy with fixed A = 2.9 GeV, which indicates that the solutions are sensitive to coupling constants. Meanwhile one 
finds that one sigma exchange interaction induces significant effects. There exist solutions with the binding energies 
around — 10 ~ MeV with appropriate coupling constants and cutoff, which are also shown in Table HVl We present 
the radial wave function R{r) and function x(r) = rR{r) for the Di — D* system with J = in Fig. [5] (b) and Fig. [7] 
(b), (c) respectively. 

One finds that big g ■ g" , small g^ ■ g'l and big h'^ arc beneficial to large binding energy. The results indicate that 
the larger the cutoff A is, the deeper the binding. 



VI. BOTTOM ANALOG 



The calculation can be easily extended to study the bottom analog of Z+(4430). For such a system, its flavor 
wavefunction is 



1 r 
1 r 



b'+b*^) + \b*+b'^ 

B+B*^) + \B*+B^) 



(26) 



(27) 



The masses of the bottom mesons are ms- = 5325.0 MeV, tub, ^ 5725.3 MeV and m^j = 5732 MeV [H, [3l|. 
Since m^j — ms* or — ms» is greater than the pion mass and less than a mass, the forms of the potentials are 
exactly those given in Table HI and [Til The difference between the charmed and bottom systems lies only in the meson 
masses. The relatively small kinematic term makes it easier to form molecular states in the bottom system. 

From the results for the D[D* system fTable HlB . one sees that the a exchange gives negligible contributions to 
the binding energy of the system. The same conclusion holds for the B'^B* system. Thus we ignore effects from a 
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Di — D* system 



J = 



9 ■ 9 


n 


n" 
9a ■ 9rT 


n-a 






J* maa; (fm) 






U.U 


U.U 


-Z.oyZ.i) ) 


i.o 


U.z 








U.2 




Z.b 


U.z 






0.58 


1.0 


-133.2(2.9) 


0.3 


0.2 


0.2 


0.55 
















n 9 




1 n 


n 9 






-0.58 


i.U 


-lOO.Ol^Z.y j 


U.o 


n 9 






0.0 


0.0 


-196.4(2.9) 


0.3 


0.2 








0.2 


-193.3(2.9) 


0.3 


0.2 






0.58 


1.0 


-4.6(2.0)7-427.0(2.9) 


1.5/0.2 


0.4/0.1 


0.6 


0.55 


















n 9 


91 7 n('9 Q"; 


n 3 


n 9 






-0.58 


i.U 


-y.i(^z.uj/-4oo.D(^^.yj 


l.i/U.z 


U.4/U.1 






0.0 


0.0 


X 


X 


X 








0.2 


X 


X 


X 






0.58 


1.0 


X 


X 


X 


-0.2 


0.55 


















0.2 












-0.58 


1.0 


X 


X 


X 






0.0 


0.0 


X 


X 


X 








0.2 


X 


X 


X 






0.58 


1.0 


X 


X 


X 


-0.6 


0.55 


















0.2 


X 


X 


X 






-0.58 


1.0 


X 


X 


X 



TABLE IV: Numerical results for the Di — D* system with several sets of parameters if we use A — 2.0 GeV and A = 2.9 GeV. 
The unit is MeV. Negative bound energies exist only for J — 0. The cross x means no bound state exists. If we use a smaller 
cutoff A = 0.7 GeV, there are no bound state solutions with the parameters in this table. 

B'l — B* system 



J = J = 1 J = 2 



9 ■ 9 


h 


E(MeV) 


rrms(fm) 


Tmax (fm) 


E(MeV) 


rrms (fm) 


rmax (fm) 


E(MeV) 


rrms(fm) 


rmax (fm) 


0.1 


0.56 


-2.6 


2.0 


1.6 


-2.6 


2.0 


1.6 


-2.8 


2.0 


1.6 


0.1 


0.84 


-14.6 


1.6 


1.5 


-14.7 


1.5 


1.5 


-14.7 


1.6 


1.5 


0.5 


0.56 


-2.3 


2.1 


1.6 


-2.5 


2.1 


1.6 


-3.2 


1.9 


1.5 


0.5 


0.84 


-14.5 


1.7 


1.5 


-14.6 


1.7 


1.5 


-14.9 


1.6 


1.5 


-0.1 


0.56 


-2.9 


2.0 


1.6 


-2.8 


2.0 


1.6 


-2.6 


2.0 


1.6 


-0.1 


0.84 


-14.8 


1.6 


1.5 


-14.7 


1.6 


1.5 


-14.7 


1.6 


1.5 


-0.5 


0.56 


-4.4 


1.7 


1.4 


-3.2 


1.9 


1.5 


-2.5 


2.1 


1.6 


-0.5 


0.84 


-15.3 


1.6 


1.5 


-14.9 


1.6 


1.5 


-14.6 


1.7 


1.5 



TABLE V: Numerical results for the B[ — B* system with several sets of parameters and A = 1.5 GeV. 



exchange in the numerical evaluation. Wc choose eight sets of parameters to calculate: \g ■ g', ft.] = [±0.1, 0.56], [±0.5, 
0.56], [±0.1, 0.84] and [±0.5, 0.84]. The results are given in Table |V] where we use the cutoff A = 1.5 GeV. With this 
cutoff, a bound state always exists with reasonable couplings. 

For the system of BiB*, we use the same sets of parameters in Table HVl since the contribution from the a exchange 
may be large. We Hst the results in Table IVll where two values of A 1.2 GeV and 1.9 GeV, are used. This case is 
similar to the case of the DiD* case. The bound states with J — may exist with the positive g ■ g" ■ 

VII. THE CUTOFF DEPENDENCE 

Up to now, we have solved the radial Schroding equation with only a few values of the cutoff A. But the binding 
energies will change with the variation of this parameter. We briefly study the cutoff dependence of the binding 
energy. 

Since a bound state is more easily formed in the bottom system, we first study the system B[B* with g ■ g' = —0.5 
and h — 0.84 for J = 0. By scanning results starting from A = 0.7 GeV, we find the binding energy has a smallest 
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Bi — B* system 



J = 



n n" 
9 9 


ll 


n n" 
9<y 9a 


n„ 




Trma (fm) 


''maa; (fm) 






U.U 


U.U 


-l.Ul^i.yj 


1 o 


n /I 
U.4 








U.2 






U.4 






0.58 


1.0 


-26.9(1.9) 


0.5 


0.3 


0.2 


0.55 
















n 9 




1.1 








-0.58 


i.U 


-oO.ol^l.y ) 


u.o 


n 9 

U.z 






0.0 


0.0 


-0.2(1.2)7-60.7(1.9) 


-iA/OA 


0.7/0.2 








0.2 


-0. 1(1. 2)/-58. 1(1.9) 


7.5/0.4 


0.7/0.2 






0.58 


1.0 


-1.3(1.2)/-111. 1(1.9) 


1.8/0.3 


0.6/0.2 


0.6 


0.55 


















n 9 


n 9^ / fi7 Q") 










-0.58 


i.U 


Ol/'l 0\/ 101 /IM o\ 

-z.n^i.zj/-izi.4(^i.yj 


1 C /A Q 


U.D/U.z 






0.0 


0.0 


X 


X 


X 








0.2 


X 


X 


X 






0.58 


1.0 


X 


X 


X 


-0.2 


0.55 


















0.2 












-0.58 


1.0 


X 


X 


X 






0.0 


0.0 


X 


X 


X 








0.2 


X 


X 


X 






0.58 


1.0 


X 


X 


X 


-0.6 


0.55 


















0.2 


X 


X 


X 






-0.58 


1.0 


X 


X 


X 



TABLE VI: Numerical results for the B\ — B* system with several sets of parameters if we use A = 1.2 GeV and A = 1.9 
GeV. Here i5(A) denotes the binding energy with corresponding cutoff. The unit is MeV. Negative bound energies exist only 
for J = 0. 

value around A — 2.3 GcV. One can always find negative eigenvalues with g ■ g' = —0.5 and h — 0.84 for this case. 
The dependence of the binding energies on the cutoff is not monotonic. We also study results with other couplings 
and J . This conclusion is always correct for the B[B* system. Some minimum binding energies with corresponding 
A are presented in Table lylTl 



Systems 


9 9 h ga ■ 9a h„ 


J=0 J=l J=2 


B[B* 


-0.5 0.84 
0.5 0.56 


-14.4(2.3) -14.0(3.5) -14.4(2.1) 
-2.3(1.4) -2.4(1.8) -2.8(2.4) 


D'lD* 


-0.5 0.84 -0.58 0.8 


-8.6(1.5) -5.9(2.1) -5.5(1.7) 



TABLE VIL Minimum binding energies E with the corresponding cutoff A for some sets of coupling constants. The unit for 
£;(A) is MeV(GeV). 

It is interesting to study whether similar behavior exists in the other systems. For the D[D* system, the minimum 
binding energy exists only for some coupling constants. One example is presented in Table IVIII For the solutions 
with h — 0.56, the binding energy decreases if we increase A until the D[D* pair is no longer bound. Thus bound 
states exist only in a small range of A. 

For the DiD* and BiB* systems, the solutions for bound states can be found only when A is big enough. The 
binding energy becomes large with the increase of the cutoff. 

In short summary, we have studied the cutoff dependence for the binding energies with some sets of coupling 
constants. Three types of behavior are found: (1) bound state solutions always exist whatever the A is. In this case, 
the binding energy reaches the minimum value at a special A; (2) bound state solutions exist with small A only. The 
binding energy becomes large when the cutoff becomes small; (3) bound state solutions exist when A is big enough. In 
this case, the binding energy increases with the cutoff. Which case is realized for a system depends on the properties 
of the components and the values of the coupling constants. One should also keep in mind the cutoff is a typical 
hadronic scale related to the system. 
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VIII. CONCLUSIONS 

QCD allows the possible existence of the glueball, hybrids, multiquarks and molecular states etc. However, none 
of them has been established firmly. Recently Belle collaboration announced the observation of charged enhancement 
Z+ (4430) , whose unique properties make Z+ (4430) hardly understood as a conventional meson. A natural explanation 
of Z+(4430) is the D[D* or DiD* molecule state. 

In this work, we re-examine the dynamics of Z"''(4430) and improve the analysis in Ref. [2] in the following aspects: 
(1) we include the sigma meson exchange contribution besides the one pion exchange potential; (2) we introduce the 
form factor to take into account the structure effect of the interaction vertex; (3) we solve the schrodinger equation 
of the S-wave D[D* or DiD* system with the help of Matlab package MATSLISE. 

We find that the one pion exchange potential from the crossed diagram plays a dominant role for the S-wave 
D[D* or DiD* system. Our numerical results indicate that with the coupling constants determined in Section |TT] 
there exists the S-wave D[D* molecular state with = 0~, 1~, 2~. The sigma meson exchange contribution to the 
binding energy is small. However one should carefully study whether the broad width of D'l disfavors the formation 
of a molecular state in the future. 

For the S-wave Di — D* system, only a molecular state with = 0~ may exist with appropriate parameters. The 
contribution from the sigma meson exchange is significant in this case. Replacing the charmed meson masses with 
bottom meson masses, we have also studied the bottom analogy Zb of the Z(4430) system. As expected, the absolute 
value of the binding energy of Zb is larger than that of Z(4430). Such a state may be searched for at Tevatron and 
LHC. Clearly future experimental measurement of the quantum numbers of Z(4430) will help uncover its underlying 
structure. 
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